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ABSTRACT 



Using the Super-Kamiokande neutrino observatory, a search was conducted 
for neutrinos produced in coincidence with gamma-ray bursts observed by the 
BATSE detector. Super-Kamiokande data in the neutrino energy range of 
7MeV ~ 100 TeV were analyzed. For gamma-ray bursts that occurred between 
1996 April and 2000 May, no statistically significant signal in excess of the back- 
ground levels was detected. Implied upper limits on associated GRB neutrino 
production are presented. 

Subject headings: GRB, Super-Kamiokande, neutrino, fluence 



1. Introduction 

Gamma-Ray Bursts (GRBs) are some of the most luminous astrophysical objects ob- 
served, and little is known about them. They have been observed to release an enormous 
amount of energy in the form of gamma-rays, and models predict that a considerable portion 
of the total energy may be carried away by neutrinos created during the burst. Unlike the 
GRB photons, which may scatter or be absorbed before they reach the Earth, the neutrinos, 
due to their small interaction cross section, arrive at the Earth virtually unaffected. Thus, 
studying neutrinos from GRBs can provide insight in understanding the processes that un- 
derlie these mysterious phenomena. This paper presents the results of a time correlation 
analysis between GRBs and Super-Kamiokande (SK) events from the neutrino data samples 
used in the solar neutrino and atmospheric neutrino analyses, as well as a direction-time 
correlation analysis using the "upward-going" muon data sample. 

SK is a 50 kton water Cherenkov detector located in the Kamioka Mine in Gifu, Japan. 
The cylindrical detector is divided into an inner and outer detector (ID and OD, respectively) 
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by a stainless-steel frame structure that serves as an optical barrier and a mounting point 
for all photo- multiplier tubes (PMTs). Cherenkov light in the ID is collected by 11,146 
inward-facing 50 cm PMTs mounted uniformly on the wall, providing 40% photo-cathode 
coverage. In the OD, 1885 outward-facing 20 cm PMTs monitor the 2.5 m thick veto region. 
The veto is used to tag incoming particles and serves as a passive shield for gamma activity 
from the surrounding rock. Detailed descriptions of the SK detector can be found elsewhere 
(Fukuda et al. 1998a; Nakahata et al. 1999). 

The most widely accepted theoretical description for gamma-ray production by GRBs 
is the relativistic fireball shock model (Rees & Meszaros 1992), in which GRBs are produced 
when relativistic ejecta from a "central engine" are slowed down by interactions, either with 
an external medium (the external shock model) or among different layers within the ejecta 
themselves (the internal shock model) (Narayan et al. 2001). The strongest confirmation 
of the fireball model comes from observations of the GRB afterglows (Katz & Piran 1997; 
Vietri 1997; Sari et al. 1998). Neutrino production in this model may be due primarily 
to p — 7 interactions (Ej, ~ 10 14 — 10 19 eV or higher) (Waxman 2000; De Paolis et al. 2001; 
Guetta et al. 2001; Waxman & Bahcall 2000; Vietri 1998) or p-n collisions (E„ - 10 MeV - 
10 GeV) (Waxman 2000; De Paolis et al. 2001; Bahcall & Meszaros 2000; Meszaros & Rees 2000; 
Kumar 1999). Great uncertainty still surrounds the "central engine" of the GRB (i.e. the 
progenitor of the fireball). The current scenarios include a collapse of massive objects (hy- 
pernovae, supernovae, supranovae, collapsars) and mergers of binary systems (black holes, 
neutron stars, white dwarfs, helium stars). Almost all models involve accretion disks and 
jet formation, which gives rise to relativistic fireballs due to extremely high temperatures 
associated with accretion. A review of these many models can be found in Meszaros (2001) 
and Piran (2000). 

Because of the lack of specificity in models, the present analysis is spectrum-independent 
and was done for all neutrino energies to which SK is sensitive (7MeV ~ 100 TeV). The low 
energy (LE) neutrino sample consists of recoil electron events from v- e elastic scattering of 
solar neutrinos, as well as background events due to radioactivity and interactions and decays 
of cosmic-ray muons. The neutrino energies in the LE sample are in the range: E u = 7 ~ 
80 MeV. More details on the LE data can be found elsewhere (Fukuda et al. 1998a). The high 
energy (HE) data sample contains fully-contained (all products of the neutrino interaction 
stop in the ID) and partially-contained (some products of the neutrino interaction penetrate 
into the OD) electron and muon events from neutrino-nucleon interactions of atmospheric 
neutrinos (Fukuda et al. 1998b). The neutrino energies in the HE sample are in the range 
E u = 200 MeV ~ 200 GeV. The upward-going muon (upmu) sample yields the highest 
energy events at SK, consisting of upward-going muons created by neutrino interactions in 
the rock beneath the detector (Fukuda et al. 1999a; Fukuda et al. 1999b). Typical parent 
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neutrino energies in the upmu sample are 10 GeV and 100 GeV for muons that stop in 
the ID and the muons that penetrate into the OD, respectively, and they span the range 
E v — 2 GeV ~ 100 TeV. Both HE and upmu sample events are predominantly due to 
atmospheric neutrino interactions. 

2. GRB-neutrino correlation analysis 

In order to search for a possible GRB neutrino signal, we conducted a time correlation 
analysis of GRBs using SK's LE and HE events, and a direction-time correlation analysis 
using SK's upmu events. The list of GRBs selected for the analysis was obtained from the 
BATSE online catalog (BATSE 2000) and the non-triggered supplement 1 to the BATSE 
catalog (Kommers et al. 2001). BATSE (Burst And Transient Source Experiment) was a 
high energy astrophysics experiment in Earth orbit on NASA's Compton Gamma-Ray Ob- 
servatory. From the official start of data taking at SK on 1996 April 1 until the end of the 
BATSE mission in 2000 May, a total of 1454 GRBs were selected to match the four years of 
operational coincidence of the two experiments. A total of 1371 GRBs were used in the LE 
and HE correlation analysis, because the LE and HE data began being collected two months 
later (1996 May 31). 

The goal of all three (LE, HE, and upmu) correlation analyses was to search for an excess 
in the number of events correlated with GRBs above expected background. All events in 
both the LE and HE samples were considered background for this analysis, since a possible 
GRB signal would provide a negligible contribution to the sample. The LE and HE mean 
background rates were assumed to be constant and were measured to be (79.8±0.3) x 10~ 5 s^ 1 
and (9.9±0.1) x 10~ 5 s -1 , respectively. In both cases no statistically significant time variation 
of the background rate was found. The background for the upmu sample was calculated by 
a Monte-Carlo simulation for each GRB separately, because there is directional variation of 
the events. 

The three analyses were naturally separated into two sections. The LE and HE time 
correlation analyses are presented together because their background rates can be calculated 
from the data, and they both use the official fiducial volume of the detector for data collection. 
Since it is not always possible to infer the direction of the parent neutrino from the direction 
of the event, only a time correlation analysis was nominally performed, and a directional 
correlation was checked for the time-correlated events. In contrast, in the upmu sample 



1 These include bursts which did not activate the real-time burst detection system onboard the BATSE 
spacecraft and were discovered during an off-line search. 
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the background rates were calculated by a simulation, and a much bigger fiducial volume 
(including the surrounding rock) was used. Also, because of a strong directional correlation 
between the parent neutrino and the resulting muon, a direction-time correlation analysis 
was possible for the upmu data sample. 

2.1. GRB Search with LE and HE Neutrinos 

Since most GRB models predict neutrinos in coincidence with the photons, and there 
is no appreciable GRB neutrino flight-time delay with respect to the photons (assuming 
m u <0.01 eV/c 2 , E u >7MeV, a ±10s window centered on the GRB time was first used 

in the analysis. This tight window would detect possible GRB neutrinos arriving virtually 
simultaneously with the onset of the gamma-ray burst. Next, a time correlation window 
of ±100 s centered on the GRB time was used because 90% of the GRBs analyzed have 
T90<100s( 2 ). A conservative ±1000 s window was also used because virtually all analyzed 
GRBs (>99.9%) have T90<1000 s. In each case the time window was taken as ± in order to 
detect neutrinos that may precede or lag the photons due to some unknown GRB process. 
Finally, a search window of 1 hr was used to scan the 24 hr period before each GRB (24 
independent searches) in order to search for neutrinos that may have been produced in the 
GRB from supernova-like (SN-like) processes, and would precede the photons. Similarly, a 
1 hr window was also used to scan the 24 hr period after each GRB in an attempt to detect 
any neutrinos that may have been created in a GRB afterglow (AG) process. 

Figure 1 shows the distribution of the number of observed signal candidate events for 
each time correlation window, compared with the expectation from Poisson fluctuations of 
the measured background for the LE (left) and the HE (right) data samples. For the SN-like 
and AG searches, only the 1 hr window with the largest deviation of the signal from the 
background is shown (out of 24 windows searched in each case). The shaded region of the 
background prediction in each bin represents the 1 a statistical error bars. As demonstrated 
by Figure 1, the observed distribution of events is in good agreement with expectations from 
the background. 

Table 1 summarizes the bin with the largest number of observed events for each time 
correlation window. It shows the number of GRBs in the analysis sample, the expected 
background for a single GRB, the maximum number of observed events, the number of GRB 
candidates with that number of events (in cases of one or two GRBs, their BATSE catalog 



2 T90 is the duration (measured by BATSE) required for integration of 90% of the burst fluence for each 
gamma-ray burst. 
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number is given), and the probability to get at least that many GRB candidates, given the 
expected background and the total number of GRBs. The top entry of each search window 
corresponds to the LE analysis and the bottom entry to the HE analysis. The last entry 
in the table represents the result from the upmu search analysis (see §2.2). The number 
of GRBs is different in each search, because the three data samples employ different live- 
time calculation methods. Moreover, since we require that the detector live-time around 
each GRB at least equals the length of the time search window, longer search windows 
yielded a smaller number of GRBs to be analyzed. In both the SN-like and AG searches 
the probabilities include a trials factor of 24 to account for the 24 search windows used. 
All the probabilities are consistent with the expectation from random fluctuations of the 
background. 

For some neutrino interactions, a correlation is expected between the parent neutrino 
direction and the direction of the interaction product (eg. lepton). In order to further 
examine whether the observed correlated events are possibly signal rather than background, 
a directional correlation check was performed for the GRBs with at least one LE or HE event 
in any of their search time windows. For each GRB-event correlated pair, the cosine of the 
angle between their directions (cos 0qrb) wa s calculated. Figure 2 shows the distributions 
of cos 0grb for each time correlation window. A possible directional correlation signal would 
be manifested as a rise of the distribution at cos 6grb = 1- All the distributions are flat 
within the background fluctuations, which demonstrates no directional correlation. 



2.2. GRB Search with Upward-Going Muons 

The expected background for each GRB was, as stated earlier, calculated by simulat- 
ing the atmospheric neutrino induced upward-going muon events with a 40-year equivalent 
Monte-Carlo. These Monte-Carlo simulations use Bartol atmospheric neutrino flux (Agrawal et al. 1996), 
the GRV94 parton distribution function (Gliick et al. 1995), and energy loss mechanisms of 
muons in the rock (Lipari & Stanev 1991). In the simulation we added suppression of 
flux due to neutrino oscillations with the oscillation parameters given by sin 2 29 = 1 and 
Am 2 = 2.5 x KT 3 eV 2 (Nakahata 2001), which are SK's best-fit — > u T:X oscillation param- 
eters for the data sample used. 

The direction of an upmu is strongly correlated with the incident neutrino direction. 
The RMS angle between a muon and its parent neutrino is 3.7°, assuming an atmospheric 
neutrino spectrum (harder GRB neutrino spectra would result in a smaller angular separation 
of the muons and GRBs). Taking into account this separation and the BATSE directional 
uncertainty for each GRB (average of ~ 5°), we chose a direction search cone of a 15° half 
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angle between the GRB and the muon. 

Because we use a 15° direction cut, in this analysis we only used the conservative ±1000 s 
time window to perform the search. We found only one upmu which was correlated in time 
(within ±1000 s) and direction (within 15°) with a GRB. Since there was only one correlation 
found, it is only presented as the last entry in Table 1 (see §2.1 for the description of the 
table). Again, the probability is consistent with the expectation from the background. 

In addition, we searched for upward-going muons that were correlated in direction with 
a GRB (<15°), and correlated in time with each other (<1000 s), all within a ±24 hr window 
of a single GRB, again in hope to detect any possible SN-like or AG neutrinos. No such 
clusters of muons were found. 



3. Upper Limits on Neutrino Fluence 

In the absence of any clear neutrino signal from the GRBs, 90% C.L. GRB neutrino 
and antineutrino fluence upper limits were calculated. In principle, a fluence limit, F[cm~ 2 ] 
is given by 

F= ^° (1) 

where iVgo is the 90% C.L. limit on the total number of neutrino interactions detected in SK, 
Nt is the number of interaction targets (electrons or nucleons), a(E u ) is the total neutrino 
cross section as a function of neutrino energy, e(E v ) is the detector efficiency as a function 
of neutrino energy, and \{E V ) is the neutrino energy spectrum normalized to unity. 

Without assuming the neutrino spectrum (there is no generally accepted GRB neutrino 
spectrum), it is impossible to directly calculate the neutrino fluence. A quantity we can 
calculate without assuming a spectrum is a fluence limit "Green's function", &(E U ), which 
is obtained by replacing \(E U ) by a delta function 8(E — E v ): 

$LE ' HE {Eu) = N T J a(E' u )e(E'ME u - E' u )dE' u ' (2) 



or in case of the upmu analysis: 



^90 



< ^ ) upmu(^) r p; i j (3) 

AN «f lit W< E ^} ^ " E D dE l 



where A is the average effective area for neutrino detection, r{E ll ) is the effective range of 
the muon (including the rock), and N a is Avogadro's number. For the upmu sample, §(E U ) 
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is calculated slightly differently, because the fiducial volume includes the rock surrounding 
the detector, which also results in much lower fluence upper limits. 

The fluence limit "Green's function" represents the fluence limit on monoenergetic 
neutrinos at different specific energies. In order to obtain the total integrated neutrino 
fluence limit, F, one needs to convolute $(E V ) with a particular neutrino spectrum, X(E V ): 



The 90% C.L. upper limit on the number of SK events per GRB, N 90 , was obtained 
by assuming a Poisson process with background (Montanet et al. 1994). Table 2 shows this 
number limit used in the calculation of $>{E V ) for all three data samples. The background 
value in the table represents the total cumulative number of expected background events 
for all the GRBs in the analysis for the given search window; the signal value is the total, 
cumulative number of events detected in that search window for all GRBs. In order to set 
a limit on neutrino fluence per GRB, a total number limit was first calculated from these 
combined values, and then divided by the number of GRBs in the search. For the LE and HE 
samples, the ±100 s window was a natural choice for calculating the upper limits, because 
for 90% of the GRBs, 100 s is the time required for integrating 90% of the gamma-ray signal. 
Since in the HE sample we can distinguish e-like and //-like events, the Ago upper limits are 
given for e-like (i/ e , v e ) and //-like z> M ) separately. In this way we gain about 35% in 
sensitivity compared to a result from combined HE data. For the upmu sample only one 
search window was used. 

The fluence limits presented here reflect the sensitivity of SK's solar neutrino, atmo- 
spheric neutrino, and upward-going muon official data samples to possible GRB neutrinos. 
The detection cross section for u e , v^, and at energies below 100 MeV (solar neutrino sam- 
ple) comes predominantly from neutrino-electron elastic scattering [y e~ — > v e~); the cross 
section for v e is dominated by inverse beta decay scattering {y e p — ► n e + ). All four types of 
neutrinos at energies above 100 MeV (atmosperic neutrino sample and upward-going muon 
sample) are detected at SK by neutrino-nucleon charged current quasi-elastic scattering 
[y N — > I N') and deep inelastic scattering (single, multiple, and coherent 7r production). 
In the upmu sample, however, only the interactions that produce ^ are considered, and 
therefore, we are limited to calculating the fluence limit only for and v^. 

The efficiency of SK at each neutrino energy was calculated by generating MC inter- 
actions of monoenergetic neutrinos and simulating the detector response to the resulting 
interaction products. The simulation was based on the SK solar neutrino analysis (for LE 
sample), atmospheric neutrino analysis (for HE sample), and upward-going muon analysis 




(4) 
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(for upmu sample). The fluence limit obtained from the HE sample was calculated only up 
to 200 GeV (40 GeV for u e ), because the understanding of the detector response is limited at 
these high energies. For the lowest energies of the HE analysis, the fluence limit for and 
Pju are sllightly worse than expected, because of low detector efficiency and a lower scattering 
cross section at the lowest energies of the sample. 

Figure 3 shows the fluence limit "Green's function", <&(E V ), for u e , v e , and obtained 
from SK's LE, HE, and upmu data samples. The calculated values of <£> at each neutrino 
energy for the three analyses are shown in Tables 3, 4, and 5, respectively. In the calculation 
we assume that all the neutrinos are emitted with the same flavor. Note that Table 3 (energies 
above 100 MeV) includes $^ and $p M coming from possible decay electron interactions from 
"stealth muons" , i.e. muons or antimuons created in the 105 — 140 MeV range (below the 
Cherenkov threshold) by quasi-elastic scattering of and v^. These limits are up to a factor 
of 10 3 lower than the limits for the 7 — 80 MeV range. This increase in sensitivity comes 
from the quasi-elastic scattering cross section that becomes dominant for and in the 
105 - 140 MeV range. 

In order to demonstrate SK's sensitivity to GRB neutrinos, we assumed an E~ 2 neutrino 
spectrum to calculate the 90% C.L. upper limit on u e , i? e , z/ M , and total fluence. Table 6 
shows this total fluence calculated for neutrinos produced in SK's three neutrino energy 
ranges: 7 MeV - 80 MeV (from LE sample), 200 MeV - 200 GeV (from HE sample), and 
2 GeV — 100 TeV (from upmu sample). Note that for and the fluence limit in the LE 
range does not include the "stealth muon" limits (105 — 140 MeV range). If we extended our 
LE neutrino energy range up to ~ 140 MeV, our fluence upper limits for and in Table 6 
would become 2.64 x 10 6 cm -2 and 3.48 x 10 6 cm -2 , respectively, which is better than the 
existing limit roughly by a factor of 10 2 . 

For comparison, we calculated a rough estimate of the expected neutrino fluence, F, from 
a cosmological GRB. In the calculation we assumed an E~ 2 neutrino spectrum, a distance to 
the GRB of z=l (D L 6.6 Gpc), the total neutrino energy of 10 53 erg (Kumar 1999) emitted 
isotropically, the fact that all emitted neutrinos have the same flavor, and that the energy 
range of emitted neutrinos coincides with the energy range of each data sample. Using 
these assumptions, we predict a neutrino fluence from a single GRB to be 1.4 cm -2 (LE), 
1.7 x 10~ 2 cm~ 2 (HE), 1.1 x 10 -3 cm~ 2 (upmu), which is roughly a factor of ~ 10 6 , 10 4 , and 
30 lower than our best upper limit in the three energy ranges, respectively. 

Comparing SK's fluence limits to previously published GRB neutrino fluence upper 
limits from the 1MB detector (Becker-Szendy et al. 1995), we find that our limits are better 
roughly by a factor of 10 2 for energies >1 GeV. Our fluence limit on upward-going muons per 
average burst (0.31 x 10~ 9 cm~ 2 ) is comparable to a similar limit by the MACRO detector 
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(Ambrosio et al. 1996). Our neutrino fluence limits are consistent with the recent predictions 
for SK of a fireball model with a collapsar progenitor (Nagataki et al. 2002). 

4. Conclusion 

A time and direction correlation analysis of Super-Kamiokande events with BATSE 
gamma-ray bursts was performed for SK's solar neutrino, atmospheric neutrino, and upward- 
going muon data samples. No signal in excess of the expected background fluctuations was 
found. A 90% C.L. fluence upper limits on neutrino and antineutrino emission from an 
average GRB were calculated at various neutrino energies detected in SK. We also calculated 
the total GRB neutrino fluence upper limits for an E~ 2 spectrum in SK's three energy ranges. 
We found that our limits are at least a factor of 30 higher than our rough estimate of GRB 
neutrino emission and are consistent with most model predictions. SK's limits are the most 
stringent fluence upper limits for GRB neutrinos at energies 7MeV — 100 TeV. 

The authors acknowledge the cooperation of the Kamioka Mining and Smelting Com- 
pany. The Super-Kamiokande detector has been built and operated from funding by the 
Japanese Ministry of Education, Culture, Sports, Science and Technology, the U.S. Depart- 
ment of Energy, and the U.S. National Science Foundation, with support for individual 
researchers from Research Corporation's Cottrell College Science Award. 
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Table 1. Summary of the maximum signal bin of the correlation analyses. 



Number of GRBs 



Window 


Sample 


GRBs 


Bgd. a 


Sig b 


with Max. Sig. c 


Prob. d 


±10s 


(LE) 


1086 


0.0160 


1 


(11 GRBs) 


0.96 




(HE) 


1115 


0.00198 


1 


(3 GRBs) 


0.38 


±100s 


(LE) 


1081 


0.160 


4 


GRB 970605 


0.026 




(HE) 


1111 


0.0198 


1 


(26 GRBs) 


0.21 


±1000 s 


(LE) 


1027 


1.60 


7 


GRB 991004D, NTB 961019.51 


0.40 




(HE) 


1056 


0.198 


2 


(19 GRBs) 


0.45 


SN-like: (-8hr,-7hr) 


(LE) 


1018 


2.87 


11 


GRB 980601 


0.99 




(HE) 


1037 


0.357 


5 


GRB 990202B 


0.59 


AG: (+4hr,+5hr) 


(LE) 


990 


2.87 


9 


NTB 970525.36 


~1.00 




(HE) 


1011 


0.357 


5 


GRB 990202A 


0.58 


±1000 sec; <15° 


(upmu) 


1454 


0.00096 


1 


GRB 991004D 


0.75 



a Expected background for the corresponding search time window for a single GRB. 
b Maximum signal (maximum number of observed events) for a single GRB. 

c Number of GRBs with maximum signal. For one or two GRBs, the BATSE current catalog number 
is given (BATSE 2000; Kommers et al. 2001). 

Probability to get at least the observed number of GRB candidates with maximum signal, given the 
expected background and the total number of GRBs analyzed. The probabilities for the SN-like and AG 
search window representatives include the trials factor for the number of search windows used (24). 
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Table 2. The 90% C.L. upper limit on number of SK events per GRB, Ni 



Sample 


GRBs 


Search Window 


Tot. Bgd. a 


Tot. Sig. b 


N 90 


LE 


1081 


±100s 


173 


177 


24.1 x 10~ 3 


HE e-like 


1111 


±100s 


12.3 


16 


9.5 x 10~ 3 


HE //-like 


1111 


±100s 


9.7 


14 


9.6 x 10~ 3 


upmu 


1454 


±1000s, <15° 


0.67 


1 


2.35 x 10~ 3 



a Total, cumulated expected background for all the GRBs. 
b Total, cumulated number of observed events for all the GRBs. 



Table 3. GRB neutrino ffuence "Green's function" upper limits (90% C.L.) obtained from 

the LE neutrino data sample. 



E u $„ e [cm 2 ] $p e [cm 2 ] $"Jcm 2 ] $p M [cm 2 ] 



7 MeV 


1.05xl0 9 


1.47xl0 8 


6.41 xlO 9 


6.41xl0 9 


10 MeV 


1.53xl0 8 


4.67xl0 6 


9.37xl0 8 


9.37xl0 8 


14 MeV 


6.79xl0 7 


1.90xl0 6 


4.11xl0 8 


4.11xl0 8 


20 MeV 


3.84xl0 7 


9.09xl0 5 


2.30xl0 8 


2.30xl0 8 


30 MeV 


2.28xl0 7 


4.50xl0 5 


1.37xl0 8 


1.37xl0 8 


50 MeV 


1.35xl0 7 


2.33x10 s 


7.94xl0 7 


7.94xl0 7 


80 MeV 


9.14xl0 6 


1.84x10 s 


5.28xl0 7 


5.28xl0 7 


105 MeV 






7.55xl0 4 


1.43x10 s 


120 MeV 






5.56xl0 4 


6.09xl0 4 


140 MeV 






2.33xl0 4 


3.17xl0 4 



a The values of <& Vli and <E>p M above 105 MeV come from "stealth 
muons" produced by neutrinos in the energy range 105 — 140 MeV 
(see §3 for details). 
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Table 4. GRB neutrino fluence "Green's function" upper limits (90% CL.) obtained from 

the HE neutrino data sample. 





$„Jcm 2 ] 


$p [cm 2 ] 


$„ M [cm 2 ] 


%[cm 2 ] 


0.2 GeV 


1.08xl0 4 


1.08xl0 4 


1.74xl0 6 


1.84xl0 6 


0.4 GeV 


1.25xl0 3 


3.27xl0 3 


1.83xl0 3 


3.84xl0 3 


1 GeV 


2.49xl0 2 


9.52xl0 2 


2.88xl0 2 


9.58xl0 2 


2 GeV 


7.28X10 1 


2.11xl0 2 


7.54X10 1 


2.17xl0 2 


4 GeV 


3.61X10 1 


8.58X10 1 


3.67X10 1 


8.02X10 1 


10 GeV 


1.91X10 1 


4.21X10 1 


1.99X10 1 


3.95X10 1 


20 GeV 


8.62x10° 


1.57X10 1 


8.85x10° 


1.63X10 1 


40 GeV 


4.90x10° 


8.92x10° 


4.76x10° 


1.02X10 1 


100 GeV 




6.98x10° 


3.13x10° 


5.31x10° 


200 GeV 




2.61x10° 


1.28x10° 


2.73x10° 



Table 5. GRB neutrino fluence "Green's function" upper limits (90% C.L.) obtained from 

the upmu data sample. 



K 


ij™ 2 ] 


$p M [cm 2 ] 


2 GeV 


3.08xl0 4 


6.98 xlO 4 


4 GeV 


4.12xl0 2 


7.96xl0 2 


10 GeV 


5.87x10° 


1.32X10 1 


20 GeV 


1.13x10° 


2.18x10° 


40 GeV 


2.63X10- 1 


5.22xl0 _1 


100 GeV 


4.62xl0~ 2 


7.53xl0~ 2 


200 GeV 


1.39xl0~ 2 


2.25xl0~ 2 


400 GeV 


4.50xl0~ 3 


6.22 xl0~ 3 


1 TeV 


1.12xl0~ 3 


1.83xl0 -3 


2 TeV 


4.15xl0~ 4 


6.77xl0~ 4 


4 TeV 


1.61xl0- 4 


2.23xl0~ 4 


10 TeV 


4.87xl0^ 5 


6.72xl0~ 5 


20 TeV 


2.05xl0~ 5 


2.41xl0~ 5 


40 TeV 


8.93xl0~ 6 


1.04xl0~ 5 


100 TeV 


3.08xl0~ 6 


3.08xl0~ 6 
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Table 6. GRB neutrino total fluence upper limits (90% CL.) for a E 2 spectrum in SK's 

three energy ranges. 



Ej, range 


F, o [cm- 2 ] 


F Po [cm- 2 ] 


F, M [cm- 2 ] 


Fp M [cm- 2 ] 


Prediction* 1 [cm 2 ] 


7 MeV - 80 MeV 
0.2 GeV - 200 GeV 
2 GeV - 100 TeV 


4.44 xlO 7 
1.66xl0 2 


9.52xl0 5 
2.97xl0 2 


2.65xl0 8 ( b ) 
1.39xl0 2 
3.83xl0~ 2 


2.65xl0 8 ( b ) 
3.00xl0 2 
4.96xl0~ 2 


1.4 

1.7xl0~ 2 
1.1 x 10~ 3 



a The prediction of a GRB neutrino fluence, assuming E 2 spectrum and the matching energy 
range, z=l distance, and 10 53 total energy emitted in neutrinos. 

b These limits become 2.64 x 10 6 cm" 2 and 3.48 x 10 6 cm -2 , respectively, if we extend the LE 
energy region to ~ 140 MeV, in order to include neutrinos resulting in "stealth muons" (see §3 for 
details). 
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Fig. 1. — Distributions of the number of GRB neutrino signal candidate events (points) with 
the expected Poisson fluctuation of the background for the ±10s (a), ±100s (b), and ±1000s (c) 
time correlation windows as well as the 1 hr window between 8 and 7 hr (d) prior to GRB time 
(supernova-like process) and the 1 hr window between 4 and 5 hr (e) after the GRB time (GRB 
afterglow) , respectively. The shaded region of the background prediction in each bin represents the 
1 a statistical error bars. For each time window, the left plot is the result for the LE sample, and 
the right plot for the HE sample. 
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Fig. 2. — Distributions of the directional correlation between GRBs and the neutrino signal can- 
didate events for the ±10 s (a), ±100 s (b), and ±1000 s (c) time correlation windows as well as 
the lhr window between 8 and 7hr (d) prior to GRB time (supernova-like process) and the lhr 
window between 4 and 5hr (e) after the GRB time (GRB afterglow), respectively. For each time 
window, the left plot is the result for the LE sample, and the right plot for the HE sample. 
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Fig. 3. — The "Green's function", <&{E V ), of 90% C.L. upper limits on GRB neutrino and an- 
tineutrino fluence per GRB obtained from the LE, HE, and upmu data samples, respectively. This 
"Green's function" represents the fluence upper limit for monoenergetic neutrinos and antineu- 
trinos at different energies. To obtain the physical fluence limit, one needs to convolute <&(E V ) 
with a neutrino energy spectrum. The fluence limit for u e ends at 40 GeV, because of SK's limited 
sensitivity to u e at the highest energies. 



